Symmetry of Multiferroicity in a Frustrated Magnet TbMn2 05 



J. Okamoto,! D. J. lluang,^'^'^Ec.-Y. Mou,^ K. S. Chao,^ H.-J. Lin,i S. Park,* S-W. Cheong* and C. T. Cheni 

^National Synchrotron Radiation Research Center, Hsinchu 30076, Taiwan 
^Department of Physics, National Tsing Hua University, Hsinchu 30013, Taiwan 
Department of Electrophysics, National Chiao-Tung University, Hsinchu 30010, Taiwan 
'''Rutgers Center for Emergent Materials and Department of Physics & Astronomy, 
Rutgers University, Pitcataway, New Jersey 08854, USA 
(Dated: February 6, 2008) 

Based on measurements of soft x-ray magnetic scattering and symmetry considerations, we demon- 
strate that the magnetoelectric effect in TbMn2 05 arises from an internal field determined by 
Sq X S-q with Sq belug the magnetization at modulation vector q, whereas the magneto-elastic 
effect in the exchange energy governs the response to external electric fields. Our results set funda- 
mental symmetry constraints on the microscopic mechanism of multiferroicity in frustrated magnets. 
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Materials which exhibit coexistence of magnetism and 
ferroelectricity with cross coupling, termed multiferroic- 
ity, are attractive because they offer the possibility for 
realizing mutual control of electric and magnetic prop- 
erties. The key phenomenon behind such mutual con- 
trol lies on the capability for the induction of magneti- 
zation by an electric field or of electric polarization by 
a magnetic field, known as the magnetoelectric (ME) ef- 
fect The ME effect is an important characteriza- 
tion of multiferroicity but has been poorly understood. 
The effect could be largely enhanced by the presence of 
internal fields. However, such enhancement requires the 
coexistence and strong coupling of magnetism and ferro- 
electricity (EE), which rarely happen in real materials. 
Recent discoveries of giant magnetoelectric couplings in 
frustrated magnets thus offer new opportunities for 
a thorough scientific understanding of multiferroicity as 
well as multiferroic applications. 

In frustrated magnets, such as RMnO^ and RMn^Oj, 
(i? Tb, Dy, and Ho) |, i, i, i, 0, i, i,M [HI, E Hi , 
the spontaneous electric polarization (P) appears in cer- 
tain antiferromagnetic (AF) phases. Unlike old exam- 
ples of multiferroics, the magnetoelectric couplings ex- 
hibited by these materials are gigantic, and the magnetic 
phases involved are complicated and commonly incom- 
mensurate with lattice. The magnetic transition temper- 
ature is higher than the ferroelectric one, suggesting that 
the ferroelectricity is induced by magnetic order. Fur- 
thermore, the inversion symmetry in the magnetic phases 
with ferroelectricity is broken [111, [l2| , implying that the 
magnetic order couples to odd orders of P. In addition, 
these magnets show anomalies in the temperature de- 
pendence of dielectric constant s. For PMnOa, the ferro- 
electric transition is accompanied by a magnetic transi- 
tion from incommensurate sinusoidal to spiral AF order 
11, l2|- Kenzelmann et al. have applied the Ginzburg- 



Landau theory to understand the multiferroic behavior 
ll|. In contrast, although Chapon et al. [13] found 



that the ferroelectricity in YMn205 results from acen- 



tric spin-density waves, little is known about the under- 
lying mechanism of multiferroicity in PMn205 because 
of their structural complexity. The exact relation and 
interplay between AF order and ferroelectricity in frus- 
trated magnets are unknown and remain controversial 
Il2,il3i,il4„il5i,d6,]. 

The presence of internal fields is the simplest origin 
of the cross coupling between magnetism and ferroelec- 
tricity. Microscopically, however, it is difficult to iden- 
tify them due to their weak effect. Historically, scatter- 
ing has been shown to be powerful for measuring accu- 
mulated microscopic changes and for revealing orderings 
and their relations directly. Neutron scattering 17|, for 
example, first convincingly proved the existence of the 
antiferromagnetic phase of MnO. Here, by resorting to 
soft x-ray magnetic scattering [11 and an analysis based 
on the phenomenological Ginzburg-Landau approach, we 
demonstrate that the magnetically induced ferroelectric- 
ity in TbMn2 05 arises from an internal field determined 
by Sq' X S-^ with Sq' being the magnetization at mod- 
ulation vector q. The results indicate that the non- 
collinearity of spins is essential for the existence of in- 
duced P, independent of commensurability and chirality. 
In contrast, the external electric fields alter the exchange 
coupling, yielding anomalies in the temperature depen- 
dence of the dielectric constant. 

We present measurements of soft x-ray magnetic scat- 
tering around the L3 (2^3/2 id) absorption edge of 
Mn to reveal the detailed coupling of ferroelectricity and 
AF order in TbMn205. Soft x-ray magnetic scatter- 
ing is a newly developed technique which is sensitive 
to the magnetic moment of transition-metal d electrons 
[isl [Til [20 . [Hj . allowing us to probe magnetic order 
with high sensitivity. The scattering amplitude is propor- 
tional to the magnetization S^, which is Sje^'^'^^ with 

Sj and fj being spin moments and position vectors, re- 
spectively. We performed scattering measurements with 
the elliptically polarized-undulator beamline at National 
Synchrotron Radiation Research Center (NSRRC) in Tai- 
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FIG. 1: (Color) Photon-energy dependence of HWHM and 
scattering intensity of TbMn205 with q = (i,0,|;) and a po- 
larization at 30 K. The intensity was obtained by fitting the q 
scan at each photon energy with a Lorentzian function and a 
hnear background. Inset: Scattering geometry with a and tt' 
polarizations for incident and scattered photons, respectively. 
With q near (^,0, \ ), the angle between the scattered x-ray 
^ 6° and that between the incident x-ray 
19°. 



q , 2 , 
and the a axis is 
and the c axis is - 



wan. The inset of Fig. 1 illustrates a schematic view of 
scattering geometry. A single crystal of TbMn2O5(100) 
with dimensions of (2 x 1 x 1) mm"^ was cut and pol- 
ished to achieve a mirror-like surface, followed by a high- 
temperature annealing to remove built-up strain during 
polishing. The modulation vector q is the momentum 
transferred from the materials and lies in the scattering 
plane defined by the a and c axes, i.e., q = {qajO^qc)- 
For photons of 639 eV, the intrinsic q resolution in terms 
of half width at half maxima (HWHM) is estimated to 
be of 0.001 A^^, including the effect of photon penetra- 
tion depth. The scattering results reveal that TbMn2 05 
exhibits an AF order below 42 K. Figure 1 shows the en- 
ergy dependence of scattering intensity with q = i^,0,\) 
and HWHM, which reflects the absorption of x-rays. To 
reduce self absorption, we used photons with an energy 
just below the Mn L3 edge in the following discussion on 
the temperature and polarization dependence; for exam- 
ple, the penetration depth (~ 1800 A) of 639-eV photons 
is much larger than the correlation length of AF order 
(- 800 A defined as the inverse of HWHM). 

Temperature-dependent measurements indicate that 
the AF order of TbMn205 occurs with modulation vec- 
tors {^ztSa, 0, \ + Sc), in which Sa and Sc characterize the 
incommensurability. The temperature dependence of qa 
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FIG. 2: (Color) Temperature-dependent q of AF order of 
TbMn2 05: (A) and (B) qa component measured with tt and 
a polarizations; (C) qc component. The vertical size of rect- 
angular symbols indicates the uncertainties of qa and qc- In- 
tensities normalized to that of the commensurate AF order 
at 30 K with a polarization are expressed by means of color. 
The temperature dependence of qc is without showing rela- 
tive intensities. Intensities of incommensurate AF order with 
q = (i±(5a,0, ^ -|- 5c) are denoted as I±, and of the com- 
mensurate ordering as Ic- The insets of (A) and (C) are, 
respectively, scattering intensities of qa and qc scans at se- 
lected temperatures. 



is plotted in Fig. 2(A) and Fig. 2(B) for tt and a polar- 
izations, respectively. As the temperature decreases, the 
incommensurate AF order of TbMn205 begins to develop 
at 42 K, in agreement with neutron results For the 
temperature between 37 K and 30 K, the incommensu- 
rate scattering intensity decreases monotonically; the qa 
of the incommensurate ordering moves toward 0.5. In 
contrast, the commensurate ordering appears between 
37 K and 24 K, coexisting with the incommensurate or- 
dering. The scattering intensities plotted in the insets of 
Fig. 2(A) and Fig. 2(C) demonstrate the coexistence of 
commensurate and incommensurate AF orderings, simi- 
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lar to the coexistence of commensurate and incommen- 
surate AF phases observed in YMn205 As shown in 
Fig. 3(A), the onset of spontaneous electric polarization 
is accompanied by the incommensurate-commensurate 
AF transition at 37 K, contrary to i?Mn03 0, [HI, [H. 
On further lowering the temperature, a commensurate- 
incommensurate transition occurs at 23 K; the commen- 
surate ordering disappears. 

The polarization dependence of x-ray scattering pro- 
vides the information about the direction of magnetic 
moments. Hannon et al. [isj have shown the magnetic 
moment along a direction Z probed in x-ray scattering 
is proportional to (e'* x e) • Z, where e' is the electric 
field of the scattered light. For an incident x-ray of a po- 
larization, the scattered x-ray from TbMn205 with e'||6 
(denoted as a' polarization) makes no contribution to the 
magnetic scattering, and the scattered x-ray with e'_L6 
(tt' polarization) is predominantly sensitive to the mag- 



netic moment along the a axis, i.e. 



because n* x a 



is ~ 6° away from the a axis. Conversely, if the incident 
x-ray has tt polarization, the scattered x-ray has either tt' 
or a' polarization; tt'* x tt is parallel to the b axis, whereas 
cr'* X TT is predominantly along the c direction. As neu- 
tron measurements indicate that the magnetic moments 
are in the ab plane 0, Q , the scattering with tt polariza- 
tion is predominantly sensitive to the magnetic moment 
along the 6 axis, i.e., \Sg\. 

The knowledge of 15^1 and 15^1 enables one to inves- 
tigate how P is induced by mag netization based on the 
Ginzburg-Landau approach [ll|, [22, l23|, |24|. As hinted 
from the broken inversion symmetry, there must be odd 
orders of P coupling to S'g-. Clearly, the lowest order cou- 
pling is that an internal field Ein couples to P, and the 
free energy F can be written as F = P'^ /2xo — Ein ■ P 
with xo being the electric susceptibility. The minimiza- 
tion of F thus leads to P = xoEin- From symmetry 
point of view, because changes sign under time re- 
versal, Ein must be quadratic in magnetization and con- 
tains at least two components in 5*^. Since both q and 
—q must be paired to make P uniform in real space, Ein 
must contain both and S-q: Under the inversion op- 



eration fi 



fj, Sq' becomes S-q: Since P changes 



sign, in order for F being invariant, Ein must change 
sign. Out of the quantities that characterize the mag- 
netic order and the underlying lattice, there are two pos- 
sible combinations for Ein'. u x {S^ x S^^) or {S^- S^q^q, 
where u is d, 5, c or their combinations 25], represent- 
ing the most important anisotropic direction in the spin- 
spin interaction. Since (Sq- ■ S'-g) = \S^\'^ is finite for 
any magnetic phase occurring below the Neel tempera- 
ture Tn, an internal field of the form (S^- S-^)q implies 
that transition temperature for non- vanishing P is identi- 
cal to Tn , which disagrees with experimental observation 
of TbMn205. On the other hand, u x (S"^ x S^^) does 
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FIG. 3: (Color online) Temperature dependence of AF order 
and ferroelectric properties of TbMn2 05: (A) Comparison of 
P and intensities '^^0-0-17^)^^^; (B) Comparison of 1-kHz e 
along the b axis and simulation of Xb with \Sqf = {la + I-w)- 
Here XoS'2(gic)_^= 0.01, Xofl2(gc) =-0.01, gi{qic) = 10 and 
gi{'lc)~0- {\SSq\^) is simulated with a Lorentzian distribu- 
tion centring at 37.5 K with a HWHM of 0.5 K. Note that 1^ 
and Itt are the average intensities for a and tt polarizations, 
respectively, i.e., la-/^ = (/+ + J_)/2 for incommensurate or- 
dering and /ct/tt = Ic for commensurate ordering. P and e 
are reproduced data of zero magnetic field from Ref. 0|. 



discern different magnetic phases. In particular, it van- 
ishes for coUinear or inversion-invariant magnetic phases. 
Therefore, we conclude that Ein = J2q ^Iq^ ^ (% ^ S-^) 
is the only candidate that is consistent with symmetry 
considerations for TbMn2 05. Here 7^- is some unknown 
ftmction to be determined by microscopic models. Then 
the induced polarization P is xo X]q ^T?^ ^ {^q ^ ^-q) 

For TbMn205, because Sq — 0, by taking u — q, P 
is along the h axis, consistent with experimental obser- 
vations. The choice of m = g is also consistent with the 
work by Mostovoy 24], where the proposed expression, 
after being averaged over space, also shows that P is 
determined by x S-^. In addition, checking the mag- 
nitude of P provides another justification. Since P cx 
J2q"fq\^^\\S^\ sin(0a - 0fe) ^ with (j)a/b being phases of 

S'P'''', the appearance of P requires non- vanishing 4>a — 4>b, 
reflecting the requirement of non-coUincar spin structure 
for inducing polarization. Since the change of q is small, 
7^ is almost temperature- independent. With further as- 
sumption oi(j)a—(t>b being roughly independent of temper- 
ature, we have P cx '^^{iTrla)^ ■ Figure 3(A) shows the 
comparison of P with the sum of the intensities (It^Ict)^ 
over commensurate and incommensurate magnetic order- 
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ings. Clearly, they follow each other closely, indicating 
the validity of the proposed Ein and the assumption on 
the temperature independence of 4>a — (j^b- These obser- 
vations also indicate that the coexisting incommensurate 
and commensurate orderings break the inversion symme- 
try. Furthermore, since P is an odd function in any com- 
ponent of Sq, strong magnetic fields can simply change 
the sign of just one component in and result in the 
observed reversal of P in direction 0] . 

For a complete understanding of the ferroelectricity, 
we further discuss dielectric responses of TbMn2 05, as 
summarized in Fig. 3(B). Two anomalies are observed: 
e exhibits a sharp maximum at the incommensurate- 
commensurate transition (37 K) and a step-like structure 
at the commensurate- incommensurate transition (23 K). 
As shown below, both anomalies reflect the magneto- 
elasticity of the exchange energy. Specifically, because 
the exchange energy J is sensitive to positions of atoms 
[l^ . applying E induces 5 J that depends on the change 
of polarization 5P . Experimental data indicate that only 
the component Eh along the b axis couples to the mag- 
netic order Q so that 6 J only depends on SPb and can 
be expressed as 6J = —gi{q)Eb5Pb -\- g2{q){5PbY /2 with 
gi and g2 being positive numbers that characterize the 
corresponding couplings. Here the first term contributes 
the potential energy of the electric dipole, and when com- 
bined with —E ■ SP, it corrects Ef,. The second term is 
the elastic energy that contributes a correction to the 
original elastic energy P^/2xo and changes xo- Includ- 
ing these corrections, the electric susceptibility becomes 



Xo(l + E,5i(9)(|5,P))/(l 



XoE,92{q){\S,\')) 



The expression of Xb enables one to explain the be- 
havior of e phenomenologically. We first note that, for 
the commensurate ordering at qc , our data indicates the 
magnitude of Sg'^ saturates. It is then plausible to as- 
sume that the exchange energy reaches maximum at com- 
mensurate ordering, i.e., dSJ/dSPb — 0, which enables 
one to eliminate Eb and obtains SJ = ~g2{q){SPb)'^ /2. 
This is equivalent to set gi{qc) = and 52(90) < from 
the beginning. Since (|S'g|^) can be written as the sum 
of spin fluctuations, {\5Sq\^)^ and square of moments, 
l^g-p. The enhancement of spin fluctuations near the 
incommensurate-commensurate transition explains the 
observed sharp maximum of £ at 37 K. In contrast, 
for the incommensurate ordering below 37 K, the magni- 
tude of S^is not saturated. In this case, since gi{qic) > 
and 52(9/0) > 0, Xb is thus more sensitive to the numer- 
ator. Figure 3(B) shows the comparison of the measured 
e along the b axis and our phenomenological simulation 
of Xb with jS'^p = (/^ + Itt)- Note that Xb differs from 
e by a constant. The resemblance of these two curves 
shows that the response to E indeed arises from magneto- 
elasticity of the exchange energy. 
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